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The approach which relates thermal Green functions to forward scattering amplitudes of on- 
shell thermal particles is applied to the calculation of the gluon self-energy, in a class of temporal 
gauges. We show to all orders that, unlike the case of covariant gauges, the exact self-energy of the 
gluon is transverse at finite temperature. The leading T^ and the sub-leading ln(T) contributions 
are obtained for temperatures T high compared with the external momentum. The logarithmic 
contributions have the same structure as the ultraviolet pole terms which occur at zero temperature. 

o' 

O ' I. INTRODUCTION 

(N, 

i^ . Thermal gauge field theories in the temporal gauge has been the subject of many investigations both in the imaginary 

^^ ' and in the real time formalisms ll|-|§l. One of the main advantages of the non-covariant temporal gauge is that it 
^^ , is physical and effectively ghost-free. At finite temperature, it may be considered a more natural choice, since the 
\^ • Lorentz invariance is already broken by the presence of the heat bath. It is also convenient for calculating the response 
CN ' of the QCD plasma to a chromo-electric field, since in this case only the gluon self-energy is needed (in this gauge, the 
chromo-elcctric field depends only linearly on the gauge field A°) ||l|,0|. Despite these advantages, explicit calculations 
^►^ ' of loop corrections to Green functions are known to be more complicated than in covariant gauges. One of the 
J^ , difficulty is associated with the more involved tensorial structure of the propagator (see Eq. (nl)). The other more 
p^ ' fundamental problem is how to deal with the extra poles at g ■ n = in the gluon propagator, where q is the loop 
f^ momentum and n = (ng, 0) (jIq > 0) is the temporal axial four- vector. 

^^ In the imaginary time formalism, the standard method of calculation of thermal Green functions consists in replacing 

^~~^ the Matsubara frequency sum by a contour integral in the complex plane nn- Using this approach, all Green 
^^ functions, at one-loop order, can be written as the sum of the vacuum part plus the a thermal part which involves the 
Bose-Einstein (boson loop) or Fermi-Dirac (fermion loop) statistical distributions. 

In the case of the temporal gauge, the extra poles at g • ?i = in the gluon propagator demand more care and make 

7* ' the standard contour integral method of calculation very involved. In order to deal with this problem, Leibbrandt 

Oh. and Staley [p| have developed a technique for performing perturbative calculations at finite temperature in the 

^5 temporal gauge (see also reference H for a different approach to this problem) . Einploying a special version of the 

Mandelstam- Leibbrandt prescription P^\Uj^ combined with the (.-function method pi| , they have obtained the complete 

.!^ 1/T-expansion for the one-loop self-energy component ngQ(fco — 0,k). Their results show that the leading and sub- 

^v> leading contributions to the self-energy can be unambiguously determined in the temporal gauge, and are in agreement 

; I with the standard contour integral method employed in reference M . 

Ci ^ The temperature dependent part of Green functions can be described in terms of forward scattering amplitudes 

of on-shell thermal particles of the thermal medium. This idea was described in reference |I2] and has been further 
elaborated in the Feynman gauge and shown to be very useful for determining the partition function in QCD at high 
temperature p3| . More recently this result has been generalized to the case of general covariant gauges p4| . This 
method has been derived using both the imaginary time formalism and the real time formalism up to two-loop order 
| [15| , [l6[ . There is also an interesting relation with the path integral approach |Q . 

One of the purposes of this work is to extend the forward scattering method to a class of temporal gauges. The 
main advantage of this method is that in the high-temperature limit, it is straightforward to obtain the full tensor 
structure of both the leading T^ and the sub-leading logarithmic contributions for temperatures T such that T ^ k, 
where k denotes the external four-momentum. For the leading T^ part we obtain the well known gauge invariant 
hard thermal loop result. The gauge dependent sub- leading logarithmic part shares with the previous calculations in 
general covariant gauges p3 the interesting property of having the same structure as the ultraviolet pole contributions 
which occur at zero temperature. Another purpose of this work is to study the transversality property of the thermal 
gluon self-energy. As is well known |0,O,Q, in general covariant gauges, the exact self-energy of the gluon is not 
transverse at finite temperature. 

In section II we derive, in a class of temporal gauges, the full tensor structure of the forward scattering amplitude 
associated with the gluon self-energy. In section III we show that, in contrast to the behavior in general covariant 
gauges, the exact thermal self-energy of the gluon is transverse to all orders. In section IV we compute the leading 



gauge-invariant T^ and the sub-leading InT contributions of the non-vanishing transverse structures. The later is 
shown to be simply related to the ultraviolet pole terms which occur at zero temperature. Finally, we present a brief 
conclusion in section V. 

II. THE FORWARD SCATTERING AMPLITUDE 

The diagrams which contribute to the gluon self-energy are shown in the Figure 1. We employ the Feynman rules 
for the three- and four gluon vertices given in reference |19|| . In the class of temporal gauges characterized by a gauge 
parameter a, the gluon propagator can be written as 



7^-" 
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(1) 



where the usual axial vector has been written as n = no(l,0) = uqu. 
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FIG. 1. One-loop diagrams which contribute to the gluon-self energy. Wavy and dashed lines denotes respectively gluons 
and ghosts. All external momenta are inwards. 

Before considering the more involved diagrams in Figures (la) and (lb), let us show that the finite temperature 
contribution to the diagram in Figure (Ic) vanishes. In the imaginary time formalism this contribution is proportional 
to the following integral H] 



■oc+e 
d'q I dqoN{qo) 

— ictt+e 



n ■ qn ■ (q + k) 



+ qo^ ~qo 



(2) 



where 



N{qo) 



e ^ — 1 



(3) 



is the Bose-Einstein statistical distribution and i^^ is a momentum independent tensor which is obtained from the 
explicit expression for gluon-ghost vertex (note that the gluon-ghost vertex is not equal to zero in the class of gauges 
considered here) . Using the partial fraction decomposition (for n ■ k different from zero) 
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n ■ qn ■ (q + k) n ■ k 
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n ■ q n ■ {q + k) 



and making the change of variable q ^ q — k in the second term on the right hand side, we can see that Eq. 
vanishes when one takes into account the property 

Niqo ± fco) - N{qa ± 2n7:iT) = Niqo); n = 0, ±1, ±2 • • • . 



(4) 



(5) 



This result holds independently of any explicit prescription for the poles at q^ = 0. 

Let us now consider the diagrams in Figures (la) and (lb). We are specifically interested in those contributions 
which can be expressed in terms of forward scattering amplitudes of on-shell thermal particles. Following the steps 
described in the appendix H, such an expression can be obtained closing the contour of the go-integration in the right 
hand side of the complex plane and considering only the poles from the propagator ([^) at go = M, |g + k\ — ^o and at 
9o = IqI, respectively for the diagrams in the Figures (la) and (lb). After shifting the momentum in the contribution 



from the second pole in the diagram of Figure (la) and using the property (||), we obtain the following expression for 
the forward scattering (FS) part of the gluon self-energy 




k,\i.,a - k, V. b 

(6) 
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m^Ifs (27r)3y 2\q\ 2^ 

where the factor 1/2 in front of the curly bracket takes into account the symmetry of the graphs in figures (la) and 
(lb). It is understood that the Lorentz and color indices of the thermal particles (horizontal lines) are contracted 
with the tensor given by the curly bracket of Eq. (^ . 

The full tensor structure generated by expression m) is rather lengthy. In order to express the result in a more 
suggestive form, it is convenient to use the following tensor decomposition [|l8| 



K I FS ^ ^"'(^tPI + UlP^^ + UcP^^ + UdP;^^), (7) 



where 
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(-, ^ z,ft, • Uft,^ri,,y — fv (k^ 



ZK • VbtZnKu fC [i^uUi/ T" f^l^Unj 



fc2|fc| 



<. = ^ > (8) 

is a set of four independent tensors which can be constructed using two Lorentz indices, the gluon four momentum 
k and a vector u which describes the four velocity of the heat bath. The tensors P^^'" are transverse with respect to 
the external four-momentum fc'^ and satisfy k^P^ — and fc*Pj^ j^ {i — 1,2,3), whereas the tensors P^ and P^ 
are non-transverse. 

Using the orthogonality of the tensors given in Eq. (||) as well as the normalization ^Pj'j/P^'' "^ — P^^uP^'^ ^ — 
-kP^uP'^" ^ = P/iiP'^'' -^ = 1, it is straightforward to obtain the exphcit expressions for 11-^, 11^, 11*^, 11^ from Eqs. 
(W and (0) (the tensor algebra was performed using the Maple VR3 version of the symbolic computer package HIP 
pO| ). The resulting expressions for n*^ and 11^ turn out to be equal to zero. The vanishing of the one- loop correction 
to n-^ is a general consequence of the finite temperature Slavnov- Taylor identities fl^ while the vanishing of n*^ is a 
property which is verified by the forward scattering part of the gluon self-energy. From the transversality property of 
P^j^ and P^ it follows that 

fc^n;:t|^^ = o (9) 



III. TRANSVERSALITY OF THE EXACT GLUON SELF-ENERGY 

In order to ascertain whether the transversality property of the thermal gluon self-energy holds beyond the approx- 
imation of forward amplitudes with on-shell particles, it is useful to study the consequences of the Becchi-Rouet-Stora 
identities ua on the structure of 11°^ at finite temperature. To derive these, we start from the effective action: 



d'xd'yJ^^'^{x)Gfix - y)v\y) + - / d^xd^yA>^^ix)tli{x - y)A'"'{y) + ■■■ (10) 



where the ellipsis denote ghost terms which are unimportant for our purpose |19]. J**" is the source of the BRS 
transformation, and r]'' is a ghost field (Although in our case the closed ghost loops do not contribute to Feynman 
diagrams |2|], the ghost fields need not necessarily decouple from the BRS identity, which involves open ghost lines). 
G"f' is a quantity related to the gauge transformation of the gluon field A'^"', being given to lowest order by the 

covariant derivative 5°'^d^ — gf°-^'^A^'^ |jl3]. The quadratic part f°^ is the sum of free kinetic energy, without the 
gauge fixing term and the self-energy: 



r;i{ko, k) = 5'^\k^k, ~ k'g^,) + u;l{ko, k) 

The relevant BRS invariance of the gluon self-energy can now be written as 
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I 



d^x 



5A'l{x) dJf^'^ix) 



0. 



(11) 



(12) 



Differentiating this equation with respect to A^^{y) and ri"-{z), and setting the sources and fields equal to zero leads 
in momentum space to the identity: 



G''^'=°(fco,fc)nt(fco,fc) = 



(13) 



Since, to lowest order, G^'*^" is proportional to fc^, it follows from Eq. ( fiA\ ) that the above identity implies the relation: 

k>^ii;l = {k'g^, - k^K)G^'^' - G^- ^"^n^t, (14) 

where G^ denotes the contributions of Feynman loops to G^ . 

Let us first consider the one loop order function G^^:^ , which may be represented as shown in figure 2a. 

q 
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FIG. 2. Examples of source-ghost diagrams. The traced/wave-line represents the external source and the crossed vertex is 
proportional to g'*''. 



Using Eq. (|l|) for gluon propagator and the fact that the gluon-ghost vertex is proportional to n°', one finds that 
^^f is given at finite temperature by 



G^l^;!^ ^ag'N5''\iT) ^ 



/^ — c>o 



9 



r 



(27r)3-« (g-n)[(q-f fc) ■ n\ 



(15) 



where N is the number of colors, the space-time dimension is 4 — e and go = 2T:liT. One can now employ the identity 
(0) and make the change of variable q — > q — k in the second term, to get 



fc^ 



+ CXD 



G';:^ = ag^N5''''- (iT) V 



(f-^q 1 

(27r)3-^ (q • n) 



(16) 



Since this expression is orthogonal to {k'^g^u — k^ky), it follows from Eq. (1#, that the exact gluon self-energy is 
transverse to one-loop order. 

In order to study the behavior of G^""^ in higher orders, let us consider, for example, the two-loops diagram shown 
in Figure 2b. This yields at finite temperature a contribution given by: 



+ 00 „ 

G^ = {ag^Nf5^''i^Tf Y. j 



(f-'p 1 

(27r)3~^ (p + k) -n 



p ■ n 



(17) 
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{2T:f-- {q-n)[{q + k)-n][{p~q) 



where po = (2T:iT)j and go — (2T:iT)l. 



Using partial fraction decompositions and making appropriate changes of variables, one finds that Eq. (^7|) can be 
written in the form: 
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(18) 
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(27r)3-^ {q ■ n) 



which is rather typical of the contributions associated with two-loop diagrams. 

Comparing the relation (Uq) with the one-loop result given by Eq. (O) , one can see that both expressions have the 
same g-dependence. But such contributions vanish, in consequence of the antisymmetry of (g • n)^^. Note that, thus 
far, we have made no assumptions about the specific form of n, so that this conclusion is also true in a class of axial 
gauges ^. 

In the specific case of temporal gauges where q ■ n = qono,it is worth noticing that the integrals decouple, in fact, 
from the sums. Then, we can set, in the sense of dimensional regularisation: 



Jd'-^q\q\-=0, 



LU>0 



(19) 



Furthermore, using the Mandelstam-Leibbrandt prescription to remove the singular terms at po = and qo — 0, the 
sum over I also gives zero because of the antisymmetry of qo = (2TTiT)l. 

The above arguments, which show that G^? and Cf^s are equal to zero, may be extended to higher order to 
establish that G^'"'' must generally vanish. It then follows from Eq. (|^ that the exact self-energy of the gluon is 
transverse at all temperatures. This important property implies that, to all orders, the structure functions He and 
n^) should be equal to zero. (We show in the appendix M that these results are confirmed by explicit calculations to 
one-loop order.) In consequence of this behavior, the non-linear relation among the structure functions pq] 



n 



D — 



n2 



which must hold at finite temperature in any linear gauge, is manifestly satisfied in the class of temporal gauges 



(20) 



IV. HIGH TEMPERATURE CONTRIBUTIONS 



The non- vanishing structures 11^ and 11^ , can be expanded as a series of high-temperature terms which come from 
the hard thermal loop region of large q, such that q ^ T >> k. (There are also other contributions, such as those 
proportional to T, which originate from the region of small q iy,|6|.) This can be done in a straightforward way 
expressing the integrands in terms of the dimensionless four-vector Q = (l,q/\q\) and performing the expansion in 
powers of l/\q\- The results up to sub- leading contributions are given by 



and 
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2fc8 4fc6 



(g • fc)4 (Q . ky 



(4fc'*fc2 -I- 4fc2(Q . fc)2fc2 _ 8Q . ;.^2^3^ 



k1 



Nim 



(21) 
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where J d^ denotes the integration over the directions of q/\q\. Using that J d|g||g|A''(|g|) = 7r^T^/6 and performing 
the angular integrations, we obtain 
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(24) 



The above contributions proportional to T^ constitute the well known gauge independent hard thermal loop result 
which has been previously obtained in the Feynman gauge as well as in general covariant gauges 12^,0 ■ It also agrees 
with the results obtained in the static and long- wavelength limits using the temporal gauge with a = ||l|,^ . 

From the second lines in Eqs. ( p3| ) and (|2j) one can extract the part proportional to InT which comes from the region 

of high internal momenta. Such contributions can be readily obtained noticing that / -rYN{\q]) = —^\n{T) + • ■ •. 

In this way, inserting the logarithmic part of Eqs. (|2^ ) and (gj) into Eq. (0), we obtain the following result for the 
logarithmic contribution to the gluon self-energy in the temporal gauge 
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ab 
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(25) 



3n 

where we have substituted u — n/uQ and tiq — v? . 

Comparing the above logarithmic contribution with the dimcnsionally regularized zero-temperature gluon self- 
energy in the general axial gauge [g4|, we find that the coefficients of ln(T) and of 1/e are identical up to a relative 
sign. 



V. CONCLUSION 



We have shown that, in a class of temporal gauges, the exact thermal self-energy of the gluon is transverse to 
all orders. (Except for the Feynman gauge, where the transversality has been verified only to one loop order, this 
property is not valid in the class of covariant gauges |l^].) The forward scattering amplitude associated with the 
thermal gluon self-energy has been obtained, to one-loop order, using the approximation which consists in neglecting 
the contributions from the prescription poles. Our explicit calculations of the T^ and ln(T) terms, together with 
the arguments in jg], show that such corrections affect only the contributions proportional to powers of l/T, which 
become negligible at high temperature. Furthermore, our approach preserves the transversality property of the exact 
thermal gluon self-energy. The leading T^ and the sub-leading ln(r) contributions are consistently obtained from the 
forward scattering amplitude of on-shcU thermal particles. The relation between the dimcnsionally regularized zero 
temperature gluon self-energy and the InT contribution, given in Eq. (Eq), is relevant to explain the cancellation of 
the ln(— fc2) terms between the zero temperature and the temperature dependent parts of the gluon self-energy (such 
a cancellation has also been verified by explicit calculations in covariant gauges Pq]). 
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APPENDIX A: 



In this appendix we illustrate the forward scattering amplitude technique by considering the following type of 
integral which arises in the calculation of the two-point function in the imaginary time formalism ||7| 



/ = / d=^<7 



-^iQO--\-5 
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2m 



N{z) 



t{T,p) 



{q ■ qf {p ■ pY 



+ z 



-^) 



(Al) 



where p = q + k, fco = 2rmTi, m = 0, ±1, ±2 • ■ •, q = (z, q), and we have factorized the denominators l/{q ■ qY from the 
propagators. In covariant gauges t(q;p) has no poles in the complex z-plane, while for temporal gauges Eq. ^ yields 
extra prescription poles. For the ghost loop diagram in Fig Ic we have (in the temporal gauge) the trivial situation 
i = j = 0. The two special cases i — 1, j — I and i — 1, j — corresponds to Fig la and Figs lb, respectively. In 
what follows we will consider the case j 7^ so that there will be contributions from the poles at q-q = and p-p = 0. 
The case j = can be easily derived using the same basic steps described below. 

Since the integration in (Al) is over all values of q, it is convenient to make the change of variables q ^^ — g in all 
the terms (z ^-> — z) so that 
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Factorizing the denominators in (A2) we can write 
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The z integration can be readily performed using the Cauchy theorem and closing the contour in the right hand side 
plane where there are poles at z = |(7| and z = |p| — /cq (fco is a pure imaginary quantity at this stage of the calculation). 
In this way we obtain 
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where the subscript FS is to remind us that we have only considered the mass shell poles. Performing the change of 



variables q^q — kin the second term of (A4) we can write 
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Finally, using the property N(qo ± ko) = N{qo) and taking into account the contributions q ^^ —q we obtain 
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Using p ^ q + k, the special case when i ^ j = 1 can be written as 

t{q;p) , t{-p;-q) 
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The expression inside the bracket is a typical contribution to the on-shell forward scattering amplitude as represented, 
for example, by the first two terms in the Eq. (P). 



APPENDIX B: 



Here we show explicitly that the one- loop contributions to He and II^i , in the class of temporal gauges characterized 
by the gauge parameter a, are identical to zero. We start from the general expression 
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where the integrands Ic.oik, q) are obtained projecting the contributions from the diagrams in Figs, la and lb onto 
the tensors P^v^- Using partial fraction decompositions, as in Eq. (0), we perform shifts g ^ g — fc in all the terms 
involving (q + fc) ■ u in the denominator, which are justified in the dimensional regularization scheme. After using also 
the property coth(((7o ± ko)/2T) = coth(go/2T ± iirn) ~ coth((7o/2T) as well as the invariance of the g-integral under 
q — > q± k, the rather involved expressions of Ic.oik, q) simphfy to 



Ic{k,q) = 



]qHq + ky 



8q ■ uk ■ q — k k ■ u ~ Qk ■ qk ■ u — 8(fc • q) —-: — h 4fc g 



and 



lD{k,q) = 



q'^iq + fc)2 



2k-q-iq^ + 8^J^§^-k' 
k^ 



(B2) 



(B3) 



It is interesting to note that the terms with denominators involving factors like q-u or (q+k) -u cancel after performing 
shifts, so that we do not need to worry about prescription poles. Hence the only poles which may contribute are the 
ones at go = \<l\ and q^ = \q + k\ — fco. The corresponding residues of Icoik, q) are 
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Performing the shift g — > g — fc in the Eqs. (B5) and (|B7| ) one can see that the contributions from the distinct poles 
in Icoik, g) exactly cancel each other. 
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